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ABSTRACT 


The  ability  to  directly  oversample  and  digitize  microwave  range  signals  at  an 
antenna  is  not  possible  with  current  electronic  technologies.  The  objective  for  this  thesis 
was  to  design  and  computer  model  an  optical  sampling  and  digitization  process  using  a 
mode-locked  laser  and  fiber  lattice  accumulators.  A  novel  fiber  lattice  accumulator 
design  for  integrated  optical  sigma-delta  (LA)  digital  antenna  technology  is  presented. 

The  fiber  lattice  design  uses  phase  modulation  to  produce  the  proper  interference  between 
input  and  recirculated/delayed  optical  pulses  in  order  that  they  may  coherently  combine. 
In  this  manner,  accumulation  within  the  fiber  lattice  takes  into  account  the  sign  of  a 
sampled  bipolar  antenna  signal.  The  fiber  lattice  performance  is  numerically  evaluated 
within  a  first-order  optical  LA  digital  antenna  phase  coherent  simulation.  The  initial 
computer  simulations  show  promising  results  using  lower  frequency  antenna  signals  to 
verify  optical  design  feasibility  and  operation.  Optical  results  closely  matched  all- 
electronic  simulations.  The  error  between  the  input  antenna  and  output  signals  is 
quantified,  and  proves  correct  device  performance.  All  results  show  the  first-order 
optical  LA  does  work  and  is  ready  for  experimental  construction.  The  significance  of  this 
device  will  be  its  usefulness  in  extending  high  resolution  LA  analog-to-digital  conversion 


into  the  microwave  signal  bands. 
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IX 


I.  INTRODUCTION 


A.  BACKGROUND 

Analog-to-digital  converters  (ADCs)  are  extremely  important  components  in  the 
modem  world  of  electronic  systems.  They  provide  the  critical  translation  of  the 
measured  analog  signals  into  the  digital  realm  so  that  filtering,  modification,  analysis, 
and  processing  can  be  performed.  A  partial  list  of  some  ADC  applications  includes 
process  control,  automatic  test  equipment,  video  signal  acquisition,  audio  recordings  for 
compact  discs,  and  interfaces  for  personal  computers.  Unfortunately,  the  process  of 
converting  an  analog  signal  is  limited  by  the  speed  and  resolution  capability  of  the 
sampling  and  filtering  device.  As  a  result,  a  great  deal  of  research  has  been  conducted  on 
high  speed,  high  resolution  ADCs.  One  device  that  has  found  widespread  popularity  for 
improving  signal  resolution  is  the  oversampling  sigma-delta  (EA)  ADC  filter  architecture. 

The  primary  advantages  of  the  oversampling  EA  ADC  is  its  increased  resolution 
capability  without  the  use  of  complicated  circuitry.  It  is  also  simple  to  integrate  into  most 
systems  due  to  analog  signal  processing  component  requirements  that  are  much  less  than 
the  overall  resolution  capability  of  the  converter.  In  it’s  most  basic  form,  it  only  requires 
a  1-bit  ADC. 

Oversampling  EA  analog-to-digital  (A/D)  conversion  is  the  combination  of  two 
processes  acting  together.  The  first  is  the  oversampling  of  the  analog  signal.  This  is  a 
technique  for  sampling  higher  than  the  Nyquist  rate  to  improve  resolution.  Sigma-Delta 
Modulation  (EAM)  is  a  noise  shaping  filter  which  can  further  improve  resolution  when 
applied  to  an  oversampled  signal.  It  uses  integration  and  feedback  in  iterative  loops  to 
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track  the  incoming  oversampled  signal,  and  calculate  and  subtract  off  the  quantization 
error. 

Oversampling  EA  A/D  conversion  is  a  well  documented  technique  which  has 
found  applications  in  many  areas.  Some  of  these  include  audio,  telecommunications, 
medical  imaging,  and  cellular  phones.  One  drawback,  which  has  limited  its  bandwidth 
application  in  the  past,  is  the  inability  to  obtain  oversampling  rates  sufficient  to  sample 
signals  in  the  higher  frequency  bandwidths.  To  achieve  a  given  precision,  the  input 
sample  rate  of  the  EA  ADC  must  be  significantly  higher  than  the  conversion  rate. 
Because  of  this,  its  use  has  been  restricted  mostly  to  lower  frequencies  and  smaller 
bandwidths  due  to  the  samplir  "ate  limitations  of  electronic  technology.  This  thesis 
explores  new  areas  of  EA  ADC  within  the  optical  realm  in  an  attempt  to  discover  new 
ways  of  increasing  sampling  rate  in  order  to  extend  the  bandwidth  of  the  EA  ADC  into 
the  microwave  region  of  the  rf  spectrum. 

B.  PRINCIPLE  CONTRIBUTIONS 

All-electronic  oversampling  EA  ADCs  have  been  extensively  studied  and 
documented  [Ref.  1,2].  The  optical  EA,  on  the  other  hand,  has  very  few  documented 
research  efforts  [Ref.  3,4, 5, 6].  The  initial  documentation  used  in  this  thesis  came  from 
previous  thesis  work  performed  at  the  Naval  Postgraduate  school  [Ref.  6, 7]. 

This  thesis  is  divided  into  two  distinct  parts.  The  first  consists  of  the  theory  and 
modeling  of  the  all-electronic  EA  ADC.  The  second  part  covers  the  optical  EA  digital 
antenna.  Computer  simulations  were  developed  to  model  the  electronic  and  optical  EA 
architectures.  Results  from  the  simulations  were  then  compared  to  determine  the 
feasibility  of  developing  an  optical  EA  device.  The  theory  and  application  of  mode- 
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locked  lasers[Ref.  8],  phase  modulators,  directional  couplers,  [Refs.  9,10]  and  fiber  optic 
lattice  signal  processing  [Refs.  1 1, 12]  is  described  and  integrated  into  the  optical  EA 
digital  antenna  design. 

C.  THESIS  OUTLINE 

Chapter  13  covers  the  basic  all-electronic  EA  ADC  theory  and  presents  computer 
models  and  results  to  build  a  foundation  for  comparison  to  the  optical  EA  architecture 
that  has  been  developed.  Emphasis  was  placed  on  signal  demodulation  through  the 
process  of  decimation.  The  purpose  was  to  produce  outputs  that  could  be  readily 
identified  and  compared  with  signals  obtained  from  the  optical  models.  Sinusoidal  as 
well  as  ramp  signals  were  used  to  test  the  models. 

In  Chapter  HI  the  phase  coherent  optical  EA  antenna  architecture  is  presented 
along  with  a  description  of  each  optical  component  and  it’s  theory  and  application.  A 
description  of  developed  computer  models,  results,  and  comparison  analysis  to  the  all- 
electronic  EA  device  is  provided. 

Finally,  Chapter  IV  will  offer  conclusions,  recommendations,  and  suggestions  for 
•further  study  within  this  area. 
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II.  OVERSAMPLING  ALL-ELECTRONIC  SIGMA-DELTA  (EA) 
ANALOG-TO-DIGITAL  CONVERTERS  (ADCs) 

A.  INTRODUCTION 

Sigma-Delta  (XA)  modulators  and  their  application  in  analog-to-digital  (A/D) 
conversion  were  first  proposed  in  the  1 960's,  but  due  to  their  requirement  for  high 
oversampling  rates,  did  not  prove  practical  until  recent  improvements  in  technology. 
Today,  the  basic  theory  and  operation  of  XA  modulators  is  well  documented  with 
hundreds  of  research  papers  and  results.  Widespread  application  has  also  developed  in 
many  areas  of  technology  including  cellular  phones  and  audio  systems.  This  thesis  will 
not  attempt  to  fully  analyze  the  all-electronic  XA  ADC,  but  instead,  will  review  basic 
theory,  describe  computer  modeled  simulations,  and  show  results  for  comparison  to  the 
optical  XA  digital  antenna.  A  complete  description  of  the  XA  theory  and  application  can 
be  found  in  [Refs.  1,2]. 

B.  OVERSAMPLING  XA  ADC  THEORY  AND  DESCRIPTION 

1.  Oversampling  Theory 

Oversampling  is  the  sampling  of  an  analog  signal  at  a  rate  higher  than  the  Nyquist 
criteria,  fs=2fo,  where  fs  is  the  sampling  frequency  and  fo  is  the  frequency  of  the  input 
analog  signal  to  be  digitized.  Quantization  noise  is  a  byproduct  of  the  sampling  process 
and  usually  the  dominant  source  of  error  introduced  into  the  signal  bandwidth  when  the 
signal  is  digitized.  Other  types  of  noise,  including  thermal  and  jitter  noise,  are  also 
present,  but  usually  at  lower  levels. 

When  sampling  occurs  at  the  Nyquist  rate,  the  quantization  noise  is  at  the  highest 
level.  As  the  sampling  rate  increases,  a  more  precise  representation  of  the  analog  signal 
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is  created  and  the  quantization  noise  is  reduced.  This  reduction  is  observed  as  a  lowering 
of  the  signal  noise-floor,  or  increase  in  the  signal-to-noise  ratio  (SNR)  in  the  frequency 
spectrum  [Ref.  1,2].  During  oversampling,  the  same  amount  of  quantization  noise  found 
at  the  Nyquist  rate  is  physically  spread  over  a  larger  bandwidth  associated  with  the 
oversampling  frequency  (fs).  This  forces  a  larger  portion  of  the  quantization  error  outside 
the  signal  band.  A  filter  can  now  be  used  to  remove  the  out-of-band  frequency 
components  leaving  the  signal  band  intact.  The  result  is  a  signal  band  with  an  improved 
SNR.  An  example  of  oversampling  is  shown  graphically  in  Figure  1. 


-  l/2fs  fs  2fs  3fs  4fs 

qns  (a) 


l/2fs  fs  2fs  3fs  4fs 


(C) 

Figure  1 :  In-band  quantization  noise  power  density  with  (a) 

Nyquist  sampling,  (b)  oversampling  factor  of  4,  (c) 
output  after  filtering. 

For  calculation  and  analysis  purposes,  the  amount  of  oversampling  is  usually 
quantified  as  an  oversampling  ratio  (OSR)  that  is  written 
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When  the  OSR  doubles,  the  SNR  improves  by  3  dB,  increasing  the  resolution  of  the 
device  by  1.5  bits  [Refs.  1,2]. 

Oversampling  does  have  limitations.  In  general,  by  oversampling  a  signal, 
bandwidth  is  traded  for  signal  resolution.  As  the  sampling  rate  increases  the  resolution 
increases.  Concurrently,  the  bandwidth  of  the  system  decreases. 

2.  Sigma-Delta  (EA)  Modulation 

Sigma-delta  (EA)  modulation  is  a  noise  shaping  technique  which  combines 
integration  and  feedback  in  iterative  loops  to  obtain  high  resolution  A/D  conversion.  The 
integration  consists  of  accumulating  the  quantization  error  between  the  sampled  input  and 
the  feedback  signal.  When  the  accumulated  error  reaches  a  specific  threshold  voltage, 
the  output  signal  is  quantized,  fed-back,  and  subtracted  from  the  input.  This  effectively 
moves  the  noise  component  of  the  signal  up  in  frequency  and  outside  the  signal  band.  A 
filter  is  then  used  to  remove  the  out-of-band  frequency  components  as  in  the 
oversampling  technique  discussed  in  Section  Bl. 

The  integrator  consists  of  a  sample  delay  (1  /fs),  a  feed  back  loop  which 
accumulates  the  signal,  and  two  gain  amplifiers  which  control  accumulation  rate. 
Depending  on  circuit  order  and  architecture,  the  integrator  is  modified  so  the  sample 
delay  is  either  in  the  feed-forward  or  feedback  path.  The  reasons  and  uses  for  the 
different  integrators  will  become  apparent  when  the  entire  EA  A/D  converter  is 
introduced  in  Section  B3.  Equations  2  and  3  are  the  transfer  functions  for  the  two 
separate  integrator  architectures  [Ref.  6] 

Hfb(z)  =  — (2) 

1  -Bz~l 


l 


(3) 


Hff(z)  = 


1  -  Dz~l 


where  Hfb  (z)  (Eqn.  2)  is  the  integrator  with  the  delay  function  in  the  feedback  path  and 
Hff(z)  (Eqn.  3)  is  the  integrator  with  the  delay  in  the  feed-forward  path.  Figure  2  shows 
both  integrator  block  diagrams.  Both  work  in  the  same  manner  with  their  accumulation 


(b) 


Figure  2:  Block  Diagrams  for  ZA  Modulator  Integrators  with  the  delay 

Function  in  the  (a)  Feed  Forward  Position  ( Hff(z ))  (b)  Feedback 
Position  {HFb(z)). 

rates  being  the  difference  between  them.  The  Hff  (z)  accumulates  faster  than  the  Hfb  (z). 
The  integrator,  because  it  ‘accumulates’  the  signal  error,  is  normally  referred  to  as  an 
accumulator  when  discussing  the  sampled-data  ZA  ADC. 

3.  Single-Bit  All-Electronic  Oversampling  ZA  ADCs 
The  oversampling  all-electronic  ZA  ADC  uses  oversampling,  accumulation  and 
decimation  (lowpass  filtering  and  re-sampling  at  the  Nyquist  rate)  to  reduce  the  noise 
within  the  signal  band.  It  is  essentially  a  sampled  data  circuit,  and  thus,  the  integration  is 
performed  via  accumulation  [Ref.  6].  The  analog  signal  is  assumed  oversampled  well 
above  the  Nyquist  frequency.  The  sampled  data  is  then  amplitude  analyzed  in  a 
comparator/quantizer,  which  is  then  fed-back  and  subtracted  from  the  input.  This 
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quantized  feedback  signal  forces  the  average  value  of  the  quantized  output  to  track  the 
oversampled  input  signal.  Any  difference  between  the  two  accumulates  in  the  integrator, 
eventually  correcting  for  the  input  quantization  error  [Ref.  6]. 

By  continually  calculating  and  subtracting  the  quantization  error  from  the 
oversampled  input  signal,  an  accurate,  time  delayed,  digital  representation  of  the  input  is 
created  that  contains  a  significant  SNR  improvement  over  a  regularly  oversampled  signal. 
The  SNR  improvement,  as  with  oversampling,  is  incrementally  quantified  with  the  OSR. 

For  every  doubling  of  the  OSR  the  noise  improves  3(21-1)  decibels  (dBs), 
providing  L  +  1/2  extra  bits  of  resolution.  L  is  the  order  of  the  EA  ADC  and  is 
determined  by  the  number  of  accumulators  and  feedback  loops.  A  block  diagram  of  a 
first  order  EA  ADC  is  shown  in  Figure  3a. 


(a) 


Figure  3 :  Block  diagrams  of  a  (a)  1st  order  and  a  (b)  2nd  order  all- 
electronic  EA  ADC. 

The  first  order  single-bit  EA  ADC  is  characterized  by  a  single  feedback  loop  and 
one  accumulator.  The  accumulator  is  the  integrator  Hff  (z)  (Eqn.  3).  The  second  order, 
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shown  in  Figure  3b  has  two  feedback  loops  and  two  accumulators.  The  first  accumulator 
is  the  integrator  HFB  (z)(Eqn.  2),  whereas,  the  second  accumulator  is  the  Hff  (z)  (Eqn.  3). 
The  slower  accumulator  is  required  in  the  front  so  the  device  will  maintain  a  stable 
monotonic  accumulation  rate. 

Each  additional  accumulator  improves  the  SNR  in  accordance  with  the  relation 
3(2L-1)  dBs.  For  example,  a  second  order  EA  ADC  has  order  L  =  2  and  therefore,  9  dBs 
or  2.5  extra  bits  of  resolution  occur  for  every  doubling  of  the  OSR.  As  shown  in  Figure  4 
[Ref.  2],  it  can  be  determined  theoretically  that  as  the  order  of  the  device  increases  the 


Figure  4:  RMS  noise  that  enters  the  signal  band  as  a 
function  of  the  oversampling  ratio  (OSR). 


SNR  and  signal  resolution  improve.  For  example,  a  device  of  order  L=4  should  produce 
21  dBs  or  3.5  extra  bits  of  resolution  for  every  doubling  of  the  OSR.  On  the  other  hand, 
higher  order  EA  architectures  tend  to  be  unstable  and  difficult  to  model  and  design.  For 
this  reason,  first-  and  second-order  EAs  were  studied  exclusively  in  this  thesis. 
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4.  Decimation  Filtering 

The  output  of  the  EA  modulator  consists  of  the  high  frequency  modulation,  out- 
of-band  noise  components,  and  other  circuit  related  noise  which  dominate  at  the  lower 
frequencies  [Ref.  1].  A  decimation  filter  is  used  to  remove  all  the  quantization/ 
modulation  noise  and  extract  the  digitized  signal.  It  consists  of  a  lowpass  filter  and  a 
resampling  mechanism  to  down-sample  the  signal  to  the  Nyquist  rate.  This  is  best 
accomplished  in  multiple  stages  of  filtering  and  resampling.  The  first  stage  removes  the 
modulation  noise  and  down-samples  the  signal  to  some  intermediate  sampling  frequency. 
A  second  lowpass  filter  is  then  used  to  attenuate  the  out-of-band  components  before  the 
signal  is  again  resampled  at  the  Nyquist  rate.  As  the  signal  propagates  through  the  filters 
and  resampling  stages,  the  word  length  increases  to  preserve  the  resolution  [Ref.  6].  A 
more  thorough  description  is  found  in  [Refs.  1,2]. 

The  decimation  filter  is  not  restricted  to.lowpass  filters.  Other  types  of  filters  can 
also  be  used  [Ref.  2],  For  example,  there  are  passband  and  stopband  decimation  filters. 
For  the  purpose  of  this  thesis,  only  lowpass  filters  were  considered. 

C.  MODELING  THE  ALL-ELECTRONIC  EA  ADC 

MATLAB  and  SIMULINK  software  was  used  to  model  the  first-  and  second- 
order  all-electronic  EA  ADCs.  The  SIMULINK  block  diagrams  are  similar  to  the  block 
diagram  in  Figure  3  and  are  shown  in  Appendix  A,  Figure  A1  .  Each  model  consists  of 
the  input  antenna  signal,  the  feedback  loop,  the  accumulator,  and  decimator  for 
demodulation. 
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1.  Accumulator  Design  and  Modeling  Considerations 

In  the  first-order  £A  computer  program  design  (Appendix  A,  Fig.  Ala),  the 
accumulator  is  modeled  by  a  sample  delay  (1  /fs)  function  block,  a  feedback  loop  with  a 
summation  block  adding  the  accumulator  feedback  into  the  input  signal,  and  two  gain 
amplifiers  blocks.  This  arrangement  simulates  the  Hff(z)  integrator  transfer  function 
(Eqn.  3).  The  gain  coefficients,  C  and  D  regulate  the  accumulation  rate  of  the  device. 
Gain  values  between  0.1  and  1.1  are  acceptable  for  the  D  gain  coefficient.  A  D  gain 
coefficient  higher  than  1.1  and  the  signal  will  overload  the  circuit  and  cause  nonlinear 
accumulation.  On  the  other  hand,  D  coefficients  that  are  too  low  will  not  allow  sufficient 
accumulation  resulting  in  signal  loss.  The  C  gain  coefficient  scales  the  input  signal  to  the 
desired  level.  It  also  effects  accumulation  rate,  but  it's  effects  are  not  as  evident  as  the 
dominant  D  gain  coefficient.  The  best  results  are  achieved  when  the  gain  coefficients  are 
all  set  to  1.0. 

In  the  second-order  £A  architecture  (Appendix  A,  Fig  Alb),  the  HFb(z)  transfer 
function  (Eqn.  2)  is  used  in  the  first  stage  accumulator  and  the  Hff(z)  (Eqn.  3)  is  used  in 
the  second  stage.  The  slower  accumulator  (Hfb(z)  )  is  used  in  the  first  stage  to  ensure 
circuit  stability.  Values  for  the  gain  coefficients  A,  B,  C,  and  D  are  again  between  0.1 
and  1.1.  Different  combinations  provide  varying  accumulation  rates  and  outputs.  Both 
delay  function  boxes  within  the  accumulators  of  the  second-order  are  also  sample  delays 
(\/fs)  as  in  the  first-order.  Again,  as  for  the  first-order,  the  best  results  are  achieved  when 
all  gain  coefficients  are  set  to  1. 
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2.  Output  Comparator 

The  relay  function  block  models  a  single  bit  comparator/quantizer,  which 
modulates  the  signal  and  the  noise,  resulting  in  the  noise  component  being  moved  up  in 
frequency  and  outside  the  signal  band.  The  rate  at  which  the  comparator  is  clocked  is  the 
oversampling  frequency  (1  !fs).  If  the  accumulator  output  voltage  exceeds  the  comparator 
threshold,  Vf= 0,  the  comparator  output  voltage,  Vfeedback  =  IV.  Otherwise,  the  comparator 
output  voltage,  Vfeedback  =  -IV.  The  modulated  signal  (Vfeedback)  then  feeds  back  and  is 
subtracted  from  the  oversampled  input  antenna  signal. 

The  EA  architectures  presented  in  this  thesis  are  single-bit  devices  and  therefore, 
signal  amplitudes  greater  than  +  1  cannot  be  accommodated.  Multi-level  devices  are 
modeled  by  adding  a  multi-bit  ADC  in  place  of  the  comparator/quantizer  [Ref.  2]. 

3.  Oversampling  Simulations 

The  sampling  clock  mechanism  for  oversampling  the  input  signals  was  performed 
by  computer  time  steps.  The  signal  was  produced  using  a  function  block  in  the 
SIMULINK  program.  The  number  of  time  steps,  representing  the  signal  sampling  rate 
(f.)  were  specified  depending  on  the  amount  of  oversampling  required.  For  example,  if  a 
1  Hz  signal  was  modeled  and  an  OSR  of  100  was  required,  then  using  Equation  1,  it  can 
be  determined  that  fs= 200  Hz.  Then  1/200  would  be  entered  into  the  program  as  the 
sampling  step  size.  Two  types  of  signals  were  primarily  tested:  a  ramped  signal  with  a 
range  of  ±  1 V,  and  a  sine  function  with  an  amplitude  of  IV. 
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4.  Decimation  Filter  Model 


The  modulated  comparator  output  signal  is  a  combination  of  the  digitized  input 
signal,  the  modulation  noise,  and  other  out-of-band  components.  A  decimation  model 
was  developed  in  SIMULINK  to  remove  the  noise  and  out-of-band  components  to 
recover  the  digitized  input  signal.  The  recovered  signal  was  then  used  for  comparison 
and  analysis  in  the  optical  EA  design  results.  Appendix  A,  Figure  A2  shows  the 
SIMULINK  computer  model  designed  to  demodulate  and  analyze  the  recovered  signal. 

It  consists  of  multiple  decimation  steps,  each  of  which  has  of  a  separate  lowpass  filter  and 
down-sampler.  For  programming  purposes,  large  down-sampling  steps  should  be 
avoided  due  to  inaccuracies  in  the  program.  This  makes  the  number  of  decimation  filters 
dependent  on  the  OSR.  For  example,  if  the  OSR  =  100  and  the/0  =  1  Hz,  then/,=200Hz. 
This  requires  a  down-sample  factor  of  1 00  to  reduce  the  sampling  rate  to  fNyquisr=2  Hz, 
which  should  be  calculated  in  three  decimation  steps  performed  in  series.  The  first  two 
steps  each  down-sample  by  a  factor  of  5  to  the  intermediate  frequency  of  8  Hz.  The  third 
decimation  step  then  down-samples  by  a  factor  of  4  to  the  Nyquist  rate  of  2  Hz..  A  good 
rule  of  thumb  is:  avoid  down-sampling  factors  greater  than  8.  The  lowpass  filter  design 
used  in  the  computer  models  consisted  of  a  15  th  order  finite  impulse  response  (FIR) 
lowpass  filter  with  an  applied  Hamming  window  for  sinusoid  signals. 

Once  decimated  the  signal  is  displayed  on  a  scope  to  view  the  time-domain.  Then 
the  Fast  Fourier  Transform  (FFT)  is  calculated  to  view  the  frequency  domain.  Appendix 
A  provides  the  general  decimation,  FFT,  and  program  parameters  to  achieve  good  results. 
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D.  RESULTS  AND  CONCLUSIONS 

The  signals  of  choice  for  all  analysis  was  a  low  frequency  sine  wave  usually 
around  1-4  Hz  and  a  ramp  signal  with  a  ±  1  V  range.  To  model  high  frequency  signals 
on  a  computer  is  difficult  due  to  the  large  amount  of  computer  time  steps  required  to 
obtain  accurate  results.  Therefore,  the  frequencies  were  scaled  down  to  manageable 
levels  to  demonstrate  the  concepts  and  correct  operation  of  the  DA  architecture  designs. 
Figures  5  and  6  show  the  ramped  and  sinusoidal  antenna  input  signals  and  outputs  at 
various  stages  for  the  all-electronic  first-order  2A  design.  All-electronic  second-order 
time  domain  results  are  very  similar  to  first  order  and  thus  are  shown  in  appendix  A, 
Figures  A3  and  A4.  The  ramped  signal  shows  a  simulation  time  of 200  samples  and  a 
sampling  rate  of  1  sample  per  second.  The  input  frequency  of  the  sine  signal  is  f0  =1  Hz. 
It’s  sampling  rate  is^  =  32  Hz  giving  an  OSR  of  16.  Figures  5a  and  6a  are  the  input 
input  signals  and  Figures  5b  and  6b  are  the  oversampled  input  signal  with  the  feedback 
voltage  subtracted  (y(t) = Vinpur  V feedback)-  Figures  5c  and  6c  are  the  output  of  the  HFB(z) 
accumulator  and  5d  and  6d  are  the  modulated  comparator  output  (V feedback)-  The  final 
.  digitized  signal  can  now  be  filtered  from  the  modulated  output  via  the  decimation  filter. 
Figure  7  shows  the  digitized  filtered  signal  output  for  both  a  ramp  and  sinusoid  input. 
Noise  was  added  to  the  input  signals  in  Figure  7  to  demonstrate  how  effectively  noise  is 
removed  in  the  EA  architecture.  The  actual  output  is  the  smooth  signal.  The  initial  step 
response  at  the  beginning  of  the  ramped  output  and  the  phase  delay  in  the  sinusoidal 
output  are  due  to  the  initial  charge-up  of  the  lowpass  filter  during  decimation. 
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In  conclusion,  modeling  of  the  all-electronic  2A  ADC  was  successful  in  providing 
a  baseline  of  results  that  could  be  compared  with  the  optical  £A  design  presented  in 
Chapter  III.  The  SIMULINK  computer  models  (Appendix  A)  for  the  first-  and  second- 
order  EA  worked  as  expected  and  produced  good  results. 


0  20  40  60  80  100  120  140  160  180  200 

(a) 


i - > - ■ - ■ - ■ - ■ - - 1 - . - . _ I 

0  20  40  60  80  100  12  140  160  180  200 


(d) 

Figure  5:  Electronic  first-order  EA  ramp  (a)  input  signal  (b)  with 
feedback  subtracted  (c)  after  integration/accumulation 
(d)  comparator  output. 
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Figure  6:  Electronic  first-order  DA  sinusoidal  (a)  input  signal  (b)  with 
feedback  subtracted  (c)  after  integration/accumulation  (d) 
comparator  output. 
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(b) 

Figure  7:  Electronic  DA  A/D  Converter  outputs  for  (a)  a  ramp 

signal  with  noise  (b)  a  sinusoid  with  noise.  The  smooth 
signals  are  the  output. 
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III.  INTEGRATED  OPTICAL  ZA  DIGITAL  ANTENNAS 


A.  OVERVIEW 

Optical  oversampling  (EA)  digital  antennas  use  fiber  lattice  accumulators,  Mach- 
Zehnder  interferometers  (MZIs),  and  high  pulse-repetition-ffequency  (PRF)  mode-locked 
lasers  to  oversample  wideband  signals  directly  at  the  antenna.  As  high  PRF  mode-locked 
laser  technologies  continue  to  improve  (currently  PRFs  of  100+  GHz  with  time  jitter 
<100  fs),  the  oversampling  of  wideband  signals  becomes  increasingly  feasible.  This 
gives  an  advantage  to  optical  ZA  digital  antennas  over  contemporary  all-electronic 
receivers.  These  antennas  provide  the  ability  to  digitize  wideband  signals  directly  with 
high-resolution.  The  fiber  lattice  accumulators  embedded  within  a  feedback  loop  allow 
the  time  and  amplitude  jitter  requirements  of  the  mode-locked  laser  to  be  relaxed.  Other 
advantages  include  insensitivity  to  vibration  and  electromagnetic  interference,  separation 
of  the  sampling  mechanism  from  the  antenna  signal,  and  the  elimination  of  the  need  to 
down-convert  the  signals  before  digitization.  The  driving  force  behind  optical  ZA 
antenna  technology  is  the  large-bandwidth,  high-resolution  applications  that  are  beyond 
the  capability  of  the  all-electronic  ZA  architectures  (limited  typically  to  the  audio  range). 

The  main  focus  of  this  thesis  is  to  present  a  novel  fiber  lattice  accumulator  design 
for  integrated  optical  ZA  antenna  technology.  The  fiber  lattice  design  uses  phase 
modulation  to  produce  the  proper  interference  between  input  optical  pulses  and 
recirculated  optical  pulses  in  order  that  they  may  be  coherently  combined.  In  this 
maimer,  the  accumulation  within  the  fiber  lattice  takes  into  account  the  sign  of  the 
sampled  bipolar  antenna  signal.  The  fiber  lattice  performance  is  then  numerically 
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evaluated  within  a  first-order  optical  EA  digital  antenna  phase  coherent  simulation.  The 
errors  in  antenna  performance  as  a  function  of  the  input  signal  are  also  quantified. 

B.  OPTICAL  EA  DIGITAL  ANTENNAS 

A  block  diagram  of  a  first-order  optical  EA  digital  antenna  architecture  is  shown 
in  Figure  8.  The  corresponding  computer  program  is  shown  in  Appendix  B,  Figure  B 1 . 


Figure  8:  Integrated  optical  EA  digital  antenna. 

The  optical  EA  antenna  is  a  phase  coherent  device  that  relies  on  oversampling  the  antenna 
signal  directly  at  an  MZI  with  a  high  PRF  mode-locked  laser.  Since  light  intensity  can 
only  be  positive,  two  properly  phased  MZIs  are  used  to  encode  both  the  sampled  antenna 
signal  magnitude  and  polarity  (the  sampled  signal’s  sign  or  direction  for  fiber  lattice 
accumulation).  The  MZIs  are  also  used  to  subtract  the  output  comparator’s  feedback 
signal  (V feedbaCk)  from  the  antenna  signal  (Va„te„„a).  The  accumulation  of  the  laser  pulses  is 
accomplished  within  the  fiber  lattice  structure  using  the  polarity  information  from  the 
direction  MZI.  The  output  sample  of  the  fiber  lattice  is  detected  and  amplitude  analyzed 
with  a  high-speed  comparator.  The  decimation  filter  is  then  applied  to  the  output 
comparator  signal  (demodulation). 
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C.  OPTICAL  COMPONENTS  AND  PHASE  COHERENT  MODELING 

The  following  is  a  detailed  description  of  the  phase  coherent  model  and  individual 
components  of  the  optical  2A  antenna. 

1.  Mode-Locked  Laser 

The  mode  locked  laser  is  a  high-rate  pulsed  type  laser  capable  of  oversampling 
antenna  signals  directly  at  the  interferometers.  Each  individual  pulse  is  analogous  to  a 
clock  sample  in  an  electronic  ADC.  The  sampling  characteristics  of  a  typical  mode- 
locked  laser  include  PRFs  between  10  and  100  GHz,  pulsewiths  of  0.7  ps,  and  a  light 
wavelength  of  1550  nm. 

a.  Theory  and  Description 

In  a  typical  laser  cavity,  there  are  many  modes  operating  simultaneously. 
The  relative  phase  of  these  modes  is  often  random  and  incoherent  causing  the  light  from 
the  laser  to  fluctuate  as  the  modes  interact  constructively  and  destructively.  Mode¬ 
locking  is  achieved  when  distinct  longitudinal  modes  of  a  laser,  all  having  slightly 
different  frequencies,  are  combined  in  phase.  The  peak  amplitudes  of  these  oscillating 
modes,  together  in  phase,  periodically  combine  constructively  to  form  a  mode-locked 
pulse  [Ref.  8]. 

The  pulses  are  formed  when  sine  waves  in  multiples  of  2n  are  added  in 
phase  with  one  another  so  all  their  phases  are  zero  at  the  same  spatial  locations.  When 
added  together  they  produce  a  total  field  amplitude.  Squaring  the  amplitude  gives  the 
intensity  and  characteristic  pulsed  nature  of  a  mode-locked  laser  pulse.  A  detailed 
analysis  is  found  in  [Ref.  8]. 
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b.  Phase  Coherent  Model 


To  model  the  exact  parameters  of  a  mode-locked  laser  on  a  computer  is 
difficult  due  to  the  large  amount  of  computer  time  steps  required  to  obtain  accurate 
results.  Therefore,  proportionally  scaled  down  frequencies  and  pulse  duty  cycles, 
significantly  lower  than  the  actual  laser  operating  parameters  were  used  in  the  model  to 
demonstrate,  the  concept  and  correct  operation  of  the  optical  EA  digital  antenna.  Another 
consideration  to  insure  correct  operation  of  the  optical  EA  computer  model  is  to  choose 
PRFs  that  are  multiples  of  the  computer  sample  time.  Otherwise,  the  optical  delay  in  the 
fiber  lattice  is  unable  to  generate  the  correct  delay. 

The  phase  coherent  mode-locked  laser  model  is  designed  by  simulating 
the  electric  field  inside  each  pulse  as  a  sinusoidal  signal  Acoscoit,  where  a>i  is  the 
frequency  of  the  laser  light  and  A  is  the  amplitude  of  the  pulse.  Multiplying  by  a 
rectangular  pulse  gives  the  characteristic  pulsed  waveform.  Figure  9  shows  the  expanded 
view  of  a  simulated  laser  waveform  to  visualize  pulse  design.  In  an  actual  mode-locked 


laser,  the  percentage  of  pulse  width  to  pulse  repetition  interval  (pulse  duty  cycle)  is  1%  to 
10%  depending  on  sampling  requirements.  Scaled  pulse  duty  cycles  from  1%  to  50%, 
were  successfully  simulated.  On  the  other  hand,  the  laser  light,  because  it’s  frequency  is 
over  1000  times  greater  than  the  pulse  repetition  frequency  (PRF),  was  not  scaled 
proportionally  to  the  actual  mode-locked  laser  parameters  because  of  computer 
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limitations.  Actual  model  designs  used  PRFs  from  4-100  Hz,  and  light  frequencies  from 
20-400  Hz. 

2.  Mach-Zehnder  Interferometers 

The  MZIs  are  used  to  efficiently  couple  the  rf  antenna  signal  into  the  optical 
realm,  and  subtract  the  feedback  signal  from  the  antenna  signal.  Figure  10  shows  a 

(Bias  Voltage)  (Antenna  Voltage) 

0 


■+■  (Comparator 

Feedback) 

Figure  10:  Schematic  diagram  of  an  MZI  with  different  electrodes 
on  separate  paths. 

general  block  diagram  of  an  MZI. 

a.  Theory  and  Description 

An  MZI  consists  of  a  3dB  intensity  splitter,  two  optical  fiber  waveguides, 
electrodes,  and  an  intensity  combiner.  The  3dB  intensity  splitter  separates  the  mode- 
locked  laser  pulses  into  the  two  separate  waveguides.  The  electrodes,  which  are 
positioned  along  the  optical  waveguide  lengths,  receive  the  antenna,  feedback  and  DC 
bias  voltages  that  create  electric  fields  around  the  optical  waveguides.  Both  paths  of  the 
MZI  are  affected  separately  so  that  when  the  pulses  are  recombined  at  the  MZI  output, 
constructive  and  destructive  interference  occurs.  The  voltages  are  applied  separately  to 
each  path  with  opposite  polarity  (push-pull)  configuration.  For  example,  the  antenna 
voltage  is  applied  to  one  waveguide  and  the  feedback  voltage  is  applied  to  the  other.  The 
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applied  electric  field  modulates  the  phase  of  the  incoming  laser  pulse  by  changing  the 
index  of  refraction  of  the  optical  waveguides  as  a  function  of  the  voltage.  The  two 
waveguides  are  then  recombined  at  the  output  resulting  in  an  output  laser  pulse  amplitude 
modulated  by  the  analog  voltage  (constructive  and  destructive  interference).  If  no 
voltage  is  applied  to  the  electrodes,  the  pulses  recombine  coherently  for  the  maximum 
output. 

The  output  intensity  of  the  MZI  is  a  function  of  the  phase  difference  A<j> 
between  the  two  waveguide  paths  as  [Ref.  13] 

H, «  =  ^  | + i'cos[A*v)  +  9]  (4) 

1  in  Z*  L 

where  Iout  and  Iin  are  the  output  and  input  light  intensities,  respectively.  The  phase  angle 
9  is  a  DC  bias  term  that  is  used  to  adjust  separately  the  quadrature  point  of  the 
interferometers  to  create  both  a  magnitude  and  direction  MZI  transfer  function.  The 
voltage  dependent  phase  shift  is 


GKx 


(5) 


where  v(t)-Vantenna  -  V feedback  ne  is  the  extraordinary  index  of  refraction  of  the  optical 
guide,  r33  is  the  pertinent  electro-optic  coefficient,  ris  the  electro-optic  overlap 
parameter,  G  is  the  inter-electrode  gap,  and  XL  is  the  laser  wavelength.  In  terms  of  vK , 
the  voltage  required  to  shift  the  phase  by  n  radians 


A 


(6) 


where 
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(7) 


—  GlL 

Vk  ■  2L^r  • 

For  the  MZIs  modeled  in  this  thesis,  the  laser  pulsed  light  was  assumed  to 
be  traveling  in  the  x  direction  of  the  optical  waveguide  and  polarized  in  the  z  direction. 
Therefore,  the  optical  device  is  characterized  as  a  transverse  modulator  vice  a 
longitudinal  modulator.  This  allows  a  smaller  voltage  to  effect  the  phase  change  of  n 
radians,  making  the  device  more  efficient  [Ref.  9]. 

The  type  of  optical  medium  simulated  in  the  MZI  models  was  LiNbC>3. 
The  optical  parameters  for  LiNbCh  and  various  other  optical  materials  can  be  found  in 
[Ref.  10].  LiNbC>3  MZIs  are  perfectly  suited  for  the  high  bandwidth  applications 
(40GHz)  associated  with  the  optical  ZA  digital  antenna. 

b.  Phase  Coherent  Model 

The  intensity  transfer  function  given  in  Equation  4  is  used  to  model  the 
direction  MZI.  The  magnitude  MZI,  however,  requires  a  coherent  model  since  the  phase 
of  the  laser  pulse  must  be  preserved.  This  requires  modeling  the  MZI  waveguide  paths 
separately  (one  as  a  phase  modulator  and  the  other  left  unaffected)  to  determine  the 
modulation  effects  on  the  phase  of  each  pulse  as  a  function  of  the  antenna  voltage.  A 
block  diagram  of  the  magnitude  MZI  containing  a  phase  modulator  is  shown  in  Figure 
11. 


The  electric  field  for  the  path  of  the  MZI  containing  the  phase  modulator 

is  [Ref.  10] 


where 


A 

—cos 

2 


[a>Lt  +  S(vj\ 


(8) 
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is  the  phase  modulation  term  and  v(t)=  V^nm  -  Vfeedback  is  the  total  signal  applied  to  the 
interferometer  electrodes.  Rewriting, 


Ex  =  — (cosco^cosb  -sinco^sinS  ) 
2 


(10) 


allows  the  phase  modulation  due  to  the  antenna  and  feedback  signal  to  be  applied  to  the 
laser  pulse.  The  MZI  output  is  then 


Emt  =Ex+E2  =  —  (cosco^cosS  -sin©^sinS)+— cos©^  (11) 
2  2 


where  E2  is  the  electric  field  in  the  other  (unmodulated)  waveguide  path.  Figure  12 


Phase  coherent 
output  beam 


Figure  1 1 :  MZI  with  Electro-optic  phase. 

shows  the  intensity  transfer  function  (Eqn.  4)  for  both  the  direction  and  magnitude  MZIs 
along  with  the  detected  pulse  from  the  phase  coherent  model  (Eqn.  1 1)  for  an  applied 
ramp  function  -IV  to  IV  (pulse  detection  is  discussed  in  Section  C4).  Note  that  for 
v(t)> OV,  the  direction  MZI  detector’s  output  voltage  (corresponding  to  the  intensity  lout) 
exceeds  the  direction  comparator’s  voltage  threshold,  Vj=0.5.  The  direction 
comparator’s  output  voltage  is  used  to  drive  the  phase  modulator  within  the  fiber  lattice 
(accumulation  up  or  accumulation  down). 


1 


Figure  12:  MZI  intensity  transfer  functions  (solid  line)  and  phase 
coherent  MZI  detected  output  for  (a)  direction  MZI  and 
(b)  magnitude  MZI. 


3.  Fiber  Lattice  Structure 

a.  Theory  and  Description 

The  basic  fiber  lattice  structure  consists  of  two  directional  couplers,  an 
optical  amplifier,  and  a  delay  due  to  the  length  of  the  recirculating  optical  fiber.  These 
components  form  the  optical  equivalent  of  the  accumulator  in  the  all-electronic  2A 
architecture.  The  general  block  diagram  is  shown  in  Figure  13  [Ref.  13].  Blocks  Ao  and 


Figure  13:  General  Four  Port  Fiber  Lattice  Structure. 

A]  are  the  directional  couplers,  X/  and X2  are  the  input  ports,  Yj  and  F?  are  the  output 
ports,  G  is  the  optical  gain  block,  the  lines  represent  the  optical  fiber,  and  the  z1  block 
represents  the  optical  sample  delay  (1/PRF). 

By  choosing  different  combinations  of  directional  coupler  input  and 
output  ports  in  the  fiber  lattice  structure,  one  can  select  between  a  variety  of  known 
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transfer  functions  [Ref.  11].  Two  such  combinations  reproduce  the  accumulator  transfer 
functions  in  the  all-electronic  SA  architectures  (Eqns.  2,  3).  Selecting  the  input  port  X] 
and  output  port  Y2,  while  leaving^  and  Yj  unused,  produces  the  transfer  function 


1  -a0axGz~x 


(12) 


which  can  be  related  to  Equation  2  where  A=(l-ao)(l-ai)G  and  B=aoaiG  are  the  gain 
coefficients.  Inversely,  selecting  input  portX2  and  output  port  Yu  while  leaving  X/  and  Y2 


unused,  produces  the  transfer  function 


Hn(z)  = 


(1-u0)(1-Qi)Gz~ 
\-a0axGz~ 1 


which  can  be  related  to  Equation  3  where  C=(l-ao)(l-aj)G  and  D=aoaiG.  G  is  the 


optical  amplifier  gain,  a0  and  a;  are  the  coupling  coefficients  of  the  directional  couplers, 
and  z1  is  the  optical  sample  delay.  More  information  on  stability  and  theory  of  fiber 
lattice  structures  can  be  found  in  [Refs.  11,12].  A  program  to  convert  from  all-electronic 
gain  coefficients  (A,  B,  C,  D)  to  optical  coefficients  (ao,  aj,G)  is  provided  in  Appendix  C 
along  with  a  program  to  compare  their  outputs  with  a  specified  input. 


For  the  correct  operation  of  the  fiber  lattice  structure  within  the  ZA 


architecture,  the  directional  a  iers  must  be  matched  to  the  optical  gain  to  produce  a 
monotonic  or  linear  response  ith  a  constant  input.  Figure  14  shows  three  plots  of  the 
output  of  the  fiber  lattice  structure  with  the  same  constant  input,  but  varying  optical  gain 
coefficients  (G).  Figure  14a  shows  the  output  of  a  perfectly  tuned  fiber  lattice 
represented  by  the  monotonic  increase.  Figure  14b  shows  when  the  gain  is  set  too  low 
how  the  accumulation  rate  pet  and  settles  at  a  constant  value  (accumulation  too  slow). 
Figure  14c  shows  the  exponential  rise  associated  with  a  gain  that  is  too  high 
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(a)  (b)  (c) 

Figure  14:  Output  plots  for  fiber  lattice  structure  Hi2(z)  with  a  constant  input 
showing  (a)  monotonic  increase  (b)  gain  too  low  (c)  gain  too  high. 


(accumulation  too  fast).  A  monotonic  response  of  the  fiber  lattice  structure  is  important 
to  insure  the  steady  accumulation,  up  or  down,  of  the  input  signal.  A  table  of  matching 
values  for  ao,  aj,  and  G  are  listed  in  Appendix  C,  Table  Cl  [Ref.  7].  Both  fiber  lattice 
structures,  Hn(z)  and  H2i(z),  use  the  same  coefficients.  As  with  the  all-electronic 
accumulator,  the  difference  between  the  Hi2(z)  and  H2i(z)  are  their  accumulation  rates. 

b.  Directional  Coupler  Models 

A  directional  coupler  consists  of  two  individual  optical  fiber  waveguides 
of  the  same  optical  material  passing  through  a  separate  optical  medium.  Their  ability  to 
couple  light  energy  in  specified  ratios  makes  them  a  very  important  part  of  the  fiber 


(^DC) 


Figure  1 5 :  Schematic  representation  of  a  voltage  controlled 
directional  coupler. 
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lattice  structure  and  its  operation.  A  block  diagram  of  a  directional  coupler  is  shown  in 
Figure  15.  The  index  of  refraction  of  the  optical  waveguides  are  nj=ri3.  The  optical 
medium  has  an  index  of  refraction  n2.  The  index  of  refraction  of  the  medium  (n2)  can 
either  be  fixed  permanently  or  controlled  as  a  function  of  an  applied  DC  voltage  (n2(V)). 
The  voltage  creates  an  electric  field  that  changes  the  medium’s  index  of  refraction  ( n2 )  to 
equal  the  index  of  refraction  of «;  and  «?.  This  enables  the  light  energy  from  n  j  and  «3  to 
cross  through  the  optical  medium  and  couple  into  the  opposite  waveguide.  The  ratio  of 
energy  coupled  is  governed  by  the  size  of  the  electric  field. 

All  directional  couplers  have  two  input  ports  and  two  output  ports, 
although  only  three  of  the  four  in  each  coupler  are  used  in  the  fiber  lattice  design 
application  (Fig.  15).  In  the  Hi2(z)  fiber  lattice,  the  Ao  directional  coupler  uses  both  input 
ports  to  couple  the  light  energy  into  a  single  output.  The  second  output  port  is 
terminated.  The  Aj  directional  coupler  uses  one  input  and  two  outputs.  The  Ao  and  Aj 
coupler  models  for  the  H2j(z)  accumulator  are  reversed  in  the  Hn(z).  Figure  16  shows 
the  directional  coupler  models  used  in  the  H2I(z)  and  Hn(z)  accumulator  designs. 


p/out) 


Figure  16:  Block  diagrams  of  the  (a)  A ,  and  (b)  A0  directional  coupler 

computer  models  for  the  HI2(z)  fiber  lattice  structure.  The  Aj  * 
and  A0  are  interchanged  for  the  H21(z). 


The  coupling  coefficients  (ao,  aj)  in  the  fiber  lattice  transfer  function 
control  the  energy  coupling  amount  within  the  fiber  lattice  model.  They  calculate  the 
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scaled  amplitude  of  the  incoming  pulse  to  the  amount  that  would  in  reality  couple  to  the 
opposite  waveguide.  For  example,  in  Figure  16b,  if  ao  =0.4  in  the  Ao  coupler,  then  40  % 
of  the  light  energy  entering  Pi  (in)  will  couple  to  the  P2  waveguide.  Proportionally,  in  the 
second  waveguide  (P2),  1  -ao  must  equal  0.6,  and  therefore  60%  of  it's  entering  light 
energy  will  remain  in  the  P2  waveguide  to  add  coherently  to  the  40%  coupled  from  the  Pj 
waveguide.  In  this  particular  example,  only  the  P2  waveguide  output  is  required,  and  the 
other  is  terminated.  The  SIMULINK  program  model  of  the  directional  couplers  is  shown 
in  Appendix  B,  Figure  B4.  Different  ratios  for  a0  and  aj  can  be  chosen  with  each  set 
linked  to  a  specific  gain  (G)  in  order  for  the  fiber  lattice  device  to  operate  monotonically 
(Appendix  C,  Table  Cl). 

c.  Phase  Coherent  Model 

The  phase  coherent  model  of  the  optical  EA  accumulator  consists  of  a 
phase  modulator  positioned  at  the  input  of  the  fiber  lattice  structure.  The  schematic 
diagram  of  the  phase  coherent  Hj2(z)  fiber  lattice  is  shown  in  Figure  17.  Optical 


Delay  Optical  Gain 

Figure  17:  Schematic  diagram  of  the  fiber  lattice  structure 
with  recoupling  phase  modulator  at  the  input. 

coefficients  used  to  optimize  the  fiber  lattice  in  this  thesis  were  a<f= 0.6,  aj= 0.4  and 
G=4.16.  Only  a  first-order  optical  IA  digital  antenna  (Fig.  8)  is  presented  in  this  thesis  in 
which  only  the  Hi2(z)  fiber  lattice  structure  is  utilized. 


31 


The  phase  modulator  within  the  fiber  lattice  recouples  the  direction  and 
magnitude  MZI  outputs  so  the  correct  accumulation  of  the  signal  within  the  fiber  lattice 
occurs.  This  is  accomplished  by  phase  modulating  the  fiber  lattice  input  signal 
{magnitude  MZI  output)  depending  on  the  direction  MZI  comparator  voltage  (Fig.  8). 
That  is,  if  the  detected  intensity  voltage  from  the  direction  MZI  exceeds  the  direction 
comparator  voltage  threshold,  Vf= =0.5,  the  direction  comparator  output  voltage  selected  is 
cr=  0V.  Otherwise  the  direction  comparator  output  voltage  to  the  fiber  lattice  phase 

modulator  is  a= vn ,  the  voltage  required  to  phase  shift  the  magnitude  MZI  output  by  n. 

To  calculate  the  output  of  the  fiber  lattice  phase  modulator,  the  magnitude  MZI 
output  (Eqn.  1 1),  is  first  rewritten  as 


Eout  =£,  +E2  =Re[eJ^t+s)  +  eJ^].  (14) 

The  output  of  the  fiber  lattice  phase  modulator  (Ej)  is  then  calculated  by  multiplying 
Equation  14  by  the  direction  comparator  output  so  that 


Ea  =  Re 


^W+s)  +  eJ®i' 


V' 


.  It 

J—cr 


(15) 


where  <j—  0  or  v„ .  The  real  part  of  Equation  15,  the  output  of  the  fiber  lattice  phase 
modulator,  can  then  be  written  as 


7L  71 

cos  coLt  cos(£  4 - cr)  -  sin  coLt  sin(  J  4 - cr) 


7t  .  .  n 

cos  coLt  cos — cr  -  sin  coLt  sm — cr 


(16) 


If  <7=0,  then  no  modulation  of  the  magnitude  MZI  signal  occurs,  and  both  the  input  pulse 
and  the  recirculating  pulse  enter  the  A0  directional  coupler  (Fig.  17)  and  combine 


32 


coherently  (addition).  If  cr=v„,  then  the  input  signal  from  the  magnitude  MZI  receives  a 
phase  shift  of  ;rand  the  two  pulses  add  destructively  (subtraction). 

4.  Output  Photodetector 

When  the  light  exits  the  fiber  lattice  structure,  it  is  a  simulation  of  the  electric 
field  in  the  phase  altered  laser  pulse.  The  output  light  must  be  detected  and  converted 
into  a  voltage.  To  model  a  photodetector,  the  first  step  is  to  square  the  output  pulses  as 
they  exit  the  fiber  lattice  to  determine  their  intensity.  Next,  each  pulse  is  integrated 
individually.  The  signal’s  envelope  is  then  calculated  to  determine  the  detector  output 
voltage.  The  integration  can  be  expressed  as 

r(Af)  =  /(A0=fci*,(0*  (17) 

'l 

where  Epuise  is  the  electric  field  of  the  phase  modulated  output  pulse  from  the  fiber  lattice, 
V(At)  is  the  average  detected  voltage,  I  (At)  is  the  average  pulse  intensity,  and  At  is  the 
pulse  length.  A  block  diagram  of  the  photodetector  model  is  shown  in  Figure  18.  The 
SIMULINK  model  is  shown  in  Appendix  B,  Figure  B5. 


Intensity 

Figure  1 8 :  Block  diagram  of  photodetector  model. 
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5.  Output  Comparator 

The  detected  fiber  lattice  output  signal  modulates  the  output  comparator.  If  the 
detected  fiber  lattice  output  intensity  is  greater  than  0.5,  the  comparator  state  Vfeedback= 
IV,  otherwise  the  output  comparator  voltage  Vfeedback= -W .  The  feedback  voltage  is 
then  subtracted  from  the  antenna  signal,  v(t)=  Vantenna  -  Vfeedback,  at  both  the  direction  and 
magnitude  fyJZIs. 

6.  Decimation  Filtering 

Decimating  the  modulated  output  to  obtain  the  high-resolution  digital  signal  is 
accomplished  in  the  same  way  as  the  all-electronic  ZA  architecture  (lowpass  filtered  and 
then  resampled  at  the  Nyquist  rate).  The  only  consideration  prior  to  decimation  is  the 
simulated  sample  rate  of  the  laser  pulses  (PRF)  is  significantly  less  than  the  computer 
time  steps  required  to  efficiently  run  the  program.  For  example,  if  the  PRF  of  the  laser  is 
5  Hz  and  the  simulated  light  frequency  within  each  pulse  is  20  Hz,  the  minimum  amount 
of  computer  time  steps  required  to  generate  each  pulse  would  be  at  the  Nyquist  rate  (40 
Hz)  of  the  light  frequency.  Therefore,  after  modulation,  the  computer  time  steps  must  be 
down-sampled  to  equal  the  PRF  of  the  laser  prior  to  decimation,  otherwise,  the  computer 
calculated  decimation  function  will  assume  the  time  steps  to  be  the  PRF  and  will  not 
calculate  the  output  correctly.  Following  down-sampling,  the  normal  decimation 
(lowpass  filter  and  resample)  is  calculated  as  in  the  all-electronic  output  described  in 
Chapter  H 

D.  RESULTS  FOR  FIRST  ORDER  DESIGN 

Results  for  the  optical  ZA  digital  antenna  are  compared  with  the  all-electronic 
processor  in  Figures  19  and  20.  Figures  19a  and  20a- show  both  the  ramped  and 
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Optical  All-electronic 


Figure  19:  Comparison  plot  of  a  ramp  function.  Each  plot  shows  a  corresponding  stage  in  the 
2A  architectures.  (a)  Input  ramp  plotted  with  output.  The  lowpass  filter  during 
decimation  creates  the  step  function  at  the  beginning  of  the  output  signal  due  to 
capacitor  charging,  (b)  Output  of  accumulator/integrator,  (c)  Modulated  signal  at 
output  of  comparator. 

sinusoidal  input  antenna  signals  and  the  corresponding  output  signals.  The  initial  step 
response  at  the  beginning  of  the  ramped  output  and  the  phase  delay  in  the  sinusoidal 
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Figure  20:  Comparison  plot  of  a  ramp  function.  Each  plot  shows  a  corresponding  stage  in 
the  ZA  architectures,  (a)  Input  sinusoid  plotted  with  output  (dashes),  (b) 
Output  of  accumulator/integrator,  (c)  Modulated  signal  at  output  of 
comparator. 


output  are  due  to  the  lowpass  filter  initial  charge-up  during  decimation.  The  output  of  the 
accumulator  stage  (Figs.  19b,  20b)  and  output  comparator  stage  (Figs.  19c,  20c)  are  also 
shown  for  comparison.  The  ramp  signal  was  oversampled  with  a  laser  PRF  =  4  Hz  with 
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a  duty  cycle  of  50%.  The  sinusoid  signal  (1/4  Hz)  was  oversampled  with  the  laser  PRF  = 
20  Hz  with  a  duty  cycle  of  50%.  The  higher  the  PRFs  in  the  optical  model  and  clock 
sample  times  in  the  all-electronic  architecture,  the  faster  the  modulation  about  the 
average  antenna  signal  value  decreasing  visual  plot  clarity.  For  this  reason,  as  well  as 
increased  computer  time  steps  (as  discussed  in  Sections  Cl  and  C6)  results  using  low 
PRFs  and  oversampling  ratios  (OSRs)  were  calculated  and  are  shown  so  that  a  visual 
comparison  of  the  optical  IA  architecture  results  can  be  made  with  the  associated  all- 
electronic  design.  Figure  21  shows  a  spectral  comparison  of  the  optical  and  all-electronic 


Figure  21 :  Fast  Fourier  Transform  (FFT)  of  the  optical  and  all-electronic  EA 
outputs  from  Figure  20a. 

£A  output  results  shown  in  Figure  20a.  The  extra  frequency  peaks  clearly  seen  in  the  all- 
electronic  spectrum  are  the  tones  associated  with  the  first-order  design. 

Figure  22  shows  the  error  signal  for  the  optical  EA  digital  antenna  simulation. 

The  error  signal  is  calculated  as  the  time-domain  difference  between  the  input  antenna 
signals  and  the  output  digital  signal  for  both  the  ramp  (Fig.  22a)  and  the  sinusoid  (Fig. 
22b). 
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Figure  22:  Error  difference  between  the  input  signal  and  the 
output  signal  for  a  (a)  ramp  and  a  (b)  sinusoidal 
signal. 


E.  CONCLUSION 

Simulating  the  phase  coherent  nature  of  the  optical  EA  design,  a  working/reliable 
computer  model  was  developed.  All  simulation  results  confirm  the  correct  operation  of 
the  integrated  optical  model  in  comparison  to  the  all-electronic  design.  With  sampling 
rates  from  the  mode-locked  laser  expected  to  be  100  GHz,  oversampling  ratios  of  1000  or 
more  are  expected  for  signal  digitization  in  the  microwave  range  using  this  optical  EA 
architecture. 
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IV.  LIMITATIONS,  CONCLUSIONS,  AND  RECOMMENDATIONS 


A.  LIMITATIONS 

Although  the  computer  models  were  successful,  certain  parameters  associated 
with  the  realities  of  this  type  of  device  were  not  addressed  or  simulated.  These  include 
signal  losses  associated  with  the  optical  fiber  and  components,  undesired  crosstalk  within 
the  directional  couplers,  mode-locked  laser  jitter,  the  length  of  optical  fiber  required  to 
effect  the  proper  optical  delay  within  the  fiber  lattice  accumulator,  and  the  length  of  fiber 
required  for  the  optical  amplifier.  The  laser  jitter  can  be  added  to  the  already  developed 
program  by  programming  Gaussian  noise  to  vary  the  pulse  length  and  amplitude. 

Limitations  with  the  computer  program  also  effect  the  testing  and  results  if  not 
considered  prior  to  running  simulations.  The  computer  program  requires  the  simulated 
PRF  of  the  mode-locked  laser  to  be  a  multiple  of  the  computer  time  step.  Otherwise,  the 
optical  sample  delay  calculation  in  the  fiber  lattice  structure  becomes  quickly  inaccurate, 
resulting  in  improper  constructive  and  destructive  interference  with  a  loss  of  output 
signal.  Programming  the  actual  parameters  of  a  mode-locked  laser  is  not  very  practical 
on  a  computer  due  to  the  tremendous  amount  of  computer  time  steps  required  to  create  an 
accurate  result.  Therefore,  all  mode-locked  laser  parameters  were  scaled  down  to  values 
that  could  be  easily  handled  by  the  computer.  This  method  still  demonstrated  the  correct 
concept  of  an  optical  ZA  digital  antenna. 

B.  CONCLUSIONS 

Simulating  the  phase  coherent  nature  of  the  optical  ZA  design,  a  working/reliable 
computer  model  was  developed.  All  simulations  were  developed  with  the  software 
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MATLAB  and  SIMULINK.  The  results  confirm  the  correct  operation  of  the  integrated 
optical  model  in  comparison  to  the  all-electronic  design.  The  phase  modulator  fiber 
lattice  design  demonstrated  the  ability  to  properly  couple  and  accumulate  the  magnitude 
and  polarity  signal  information  within  the  first  order  EA  device.  Only  low  frequency 
signals  were  used  to  demonstrate  the  optical  EA  presented  in  this  thesis.  However,  the 
theory  and  concepts  clearly  demonstrate  the  optical  EA  digital  antenna  will  work  in  the 
microwave  signal  bands. 

In  conclusion,  the  optical  EA  digital  antenna  has  the  ability  to  digitize  analog 
signals.  Due  to  the  high  pulse  rates  associated  with  the  mode-locked  laser  this  could  lead 
to  direct  digitization  of  microwave  signals  at  significantly  high  oversampling  ratios 
(OSRs)  resulting  in  increased  signal  resolution  and  improved  sensor  sensitivity. 

C.  RECOMMENDATIONS 

Further  efforts  include  the  continued  optimization  of  the  first-order  integrated 
opticalmodel,  as  well  as,  expansion  to  a  second-order  model  (better  behaved  limit 
cycles).  A  block  diagram  of  the  proposed  second-order  optical  EA  digital  antenna  is 
shown  in  Appendix  D.  An  analysis  of  the  laser  jitter  effects  needs  to  be  completed  in 
order  to  verify  the  expected  relaxed  jitter  requirements  predicted  with  the  EA  architecture. 
Also,  a  signal  analysis  needs  to  be  modeled  and  conducted  on  the  oversampled  all- 
electronic  and  optical  EA  outputs  at  various  oversampling  ratios  (OSRs)  to  show  the 
increased  signal-to-noise  (SNR)  and  bit  resolution  obtained  using  the  EA  filter  theory. 
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APPENDIX  A:  ALL-ELECTRONIC  EA  SIMULINK  PROGRAM 
MODELS 


This  appendix  shows  the  first-  and  second-order  all-electronic  EA  SIMULINK  models, 
their  associated  subsystems  and  a  users  guide  to  setup  and  run  the  program.  The  second- 
order  results  in  the  time  domain  are  also  shown. 


(a) 


Figure  A1 :  SIMULINK  models  of  the  all-electronic  (a)  first-order  and  (b) 
second-order  EA  ADCs  {Files:  all-elect l.mdl,  all-elect2.mdl). 


Figure  A2:  SIMULINK  subsystem  (Fig.  Al)  model  for 
decimation  and  FFT  Analysis  of  modulated 
signals 
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Figure  A3:  Output  results  for  a  second-order  all-electronic  2A  ADC  with  a 
ramped  input,  (a)  First  stage  accumulator  output,  (b)  second 
stage  accumulator  output,  (c)  modulated  output,  (d)  input  (noisy) 
plotted  with  digitized  output. 


(a) 


(b) 


'(c) 


(d) 


Figure  A4:  Output  results  for  a  second-order  all-electronic  EA  ADC  with  a 
sinusoidal  input,  (a)  First  stage  accumulator  output,  (b)  second 
stage  accumulator  output,  (c)  modulated  output,  (d)  input  (+) 
plotted  with  digitized  output. 
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USER’S  GUIDE  FOR  ALL-ELECTRONIC  ZA  ADC  SIMULINK  MODEL 


1 .  Select  from  the  SIMULINK  library  the  type  of  antenna  signal  to  be  digitized.  Enter 
in  signal  parameters  and  connect  in  place  of  the  “Input  Signal”  block  in  Figures  Ala 
or  Alb. 

2.  Select  “Gain”  parameters.  The  ideal  for  both  models  is  to  set  all  gain  parameters  to  1 . 

3.  Set  all  sample  “Delay”  blocks  to  1  /fs,  where  fs  is  the  sampling  frequency  of  the 
system.  The  sampling  frequency  is  also  the  computer  time  step  rate  at  which  the 
overall  program  is  set  to  run. 

4.  The  “Comparator”  (Figs.  Ala,  Alb)  is  a  SIMULINK  relay  function  block.  The 
threshold  is  set  at  0.  It  should  be  set  to  output  a  “1”  when  the  block  is  “on”  and  a  “-1” 
when  the  block  is  off. 

5.  The  “Decimation”  SIMULINK  subsystem  block  (Figs.  Ala,  Alb)  lowpass  filters  and 
resamples  the  signal  to  remove  the  modulation  and  out-of-band  noise.  It  also 
calculates  the  FFT  for  the  input  and  output  signals  for  comparison.  Both  models  use 
the  same  decimation  subprogram  (Fig.  A2).  The  “FIR  decimation”  for  the  input 
signal  downsamples  the  original  oversampled  signal  to  the  Nyquist  rate.  The  “FIR 
Decimation”  block  parameters  for  the  input  and  modulated  signal  should  be  the  same. 
The  “sample  time”  is  l/fs.  The  “decimation  factor”  is  the  integer  amount  the  signal  is 
to  be  downsampled.  Avoid  downsampling  factors  greater  than  8  due  to  program 
inaccuracies.  Instead  use  multiple  decimation  blocks  in  series.  The  lowpass  filter  is  a 
“firl”.  Filter  orders  of  15-25  are  sufficient  for  sinusoidal  signals. 

6.  Following  decimation,  the  signal  can  either  be  viewed  in  the  time-domain  or  the 
frequency  domain.  In  the  frequency  domain,  a  “Hamming”  window  function  block 
and  “Buffered  FFT  Scope”  are  used  with  parameters: 

a.  Parameters  typically  used  for  the  Hamming  window  function  were:  Window 
length=128  with  periodic  sampling. 

b.  Parameters  typically  used  for  the  buffered  FFT  scope  were:  Buffer  size- 128, 
Buffer  overlap-64,  FFT  length-128,  sample  time-1#,,  where#,  equals  the 
Nyquist  frequency. 
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APPENDIX  B:  OPTICAL  2A  SIMULINK  PROGRAM  MODEL 


This  appendix  shows  the  entire  SIMULINK  optical  EA  program  design,  it’s  associated 
subsystems,  and  a  users  guide  to  setup  and  run  the  program. 


MZI  Intensity  Fnc 


Input  Signals 


ComPuter  DedSStion 
downsample 


Gain 


Figure  Bl:  SIMULINK  program  model  for  the  integrated  optical  EA  digital  antenna. 
(File:  opticall.mdT) 
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Figure  B2:  SIMULINK  subsystem  model  of  the  “Phase  Coherent  Fiber  Lattice” 
(Fig.Bl). 


_ w 

r> 

In 

Out2 

CD~* 

In 

Out 

Ini 

Coupler  A1  Coupler  AO 


Transport 

Delay 


<3* 

Gain 


Figure  B3 :  SIMULINK  subsystem  model  of 
“ Hi2(z )  Accumulator”  (Fig.  Bl) 


at 


Figure  B4:  SIMULINK  subsystem  models 
of  the  Hn(z)  directional  couplers 
(a)  Ao  and  (b)  A  /  (Fig.  B3). 
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Figure  B5:  SIMULINK  subsystem  model  the  “Photodetector” 
design  (Fig.  B2) 
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Figure  B6:  SIMULINK  subsystem  model  of  the  fiber  lattice  “Phase 
Modulator”  (Fig  B2).  Represents  Equation  16. 


3dB 

Splitter 


Figure  B7:  SIMULINK  subsystem  model  of  the  “Phase 
Coherent  Magnitude  MZI”  (Fig.  Bl). 


sin(f) 


Figure  B8:  SIMULINK  subsystem 
model  of  the  “Mode- 
Locked  Laser”  (Fig.  Bl) 
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USERS  GUIDE  FOR  OPTICAL  EA  SIMULINK  MODEL 


1 .  Select  from  SIMULINK  the  type  of  antenna  signal  to  be  digitized.  Enter  in  desired 
parameters  and  connect  in  place  of  the  “Input  Signal”  block  in  Figure  B1 

2.  Go  to  the  “Mode-Locked  Laser”  subsystem  block  (Fig.  B 1)  select  the  laser  PRF  in  the 
block  “Laser  Pulse  Window”(Fig.  B8).  Next,  choose  a  light  frequency  and  enter  into 
the  “sin(f)”  and  “cos(f)”  function  blocks. 

3 .  The  MATLAB  Function  block  labeled  “MZI  Intensity  Fnc”  (Fig.  B 1)  is  a  MATLAB 
function.  The  name  of  the  MATLAB  function  should  be  entered  in  the  SIMULINK 
block.  The  program  code  is  in  Appendix  E,  Program  E2.  The  program  should  not 
require  any  changes  unless  the  optical  properties  of  the  MZI  are  changed.  Because  it 
is  the  intensity  transfer  function,  the  only  input  required  is  the  antenna  signal  with 
feedback  already  subtracted. 

4.  The  “Phase  Coherent  Magnitude  MZI”  SIMULINK  subsystem  (Fig.  B 1)  has  the  two 
laser  outputs  as  it’s  inputs  along  with  the  antenna  signal  with  feedback  subtracted. 
Because  the  laser  sinusoidal  components  (cosine  and  sine)  are  required  to  be 
separated  in  order  to  calculate  the  second  phase  modulation  prior  to  the  fiber  lattice 
(Fig.  B7),  only  S(v)  is  calculated  and  the  rest  is  carried  over  to  the  fiber  lattice  phase 
modulator  to  perform  the  whole  calculation  at  once  before  the  signal  enters  the  fiber 
lattice  (Eqn.  16,  Chapt.  III).  Appendix  F  shows  the  SIMULINK  programs  of  the 
phase  coherent  MZI  outputs  of  both  the  magnitude  and  direction  MZIs  for 
comparison  to  the  MZI  intensity  transfer  function. 

5.  The  “Direction  Comparator”  (Fig.  Bl)  is  a  SIMULINK  relay  function  block.  The 
threshold  should  be  set  at  0.5.  It  should  be  set  to  output  “0”  when  the  block  is  “on” 
and  a  “1”  when  the  block  is  “off’. 

6.  The  “Phase  Coherent  Fiber  Lattice”  SIMULINK  subsystem  (Fig.  Bl)  consists  of  the 
“Phase  Modulator”(Fig.  B6),  “Hnfz)  Accumulator”(Fig.  B3),  and  the  “Photodetector” 
(Fig.  B5)  SIMULINK  subsystems  (Fig.  B2). 

a.  No  changes  in  the  phase  modulator  are  required  unless  a  change  in  the  optical 
parameters  is  desired.  The  value  for  v„  should  be  entered  in  the  box  labeled 
“delta”  (Fig.  B6). 

b.  The  “Gain”  in  the  accumulator  subsystems  (Fig.  B3)  must  be  matched  to  the 
directional  couplers.  The  directional  couplers  are  shown  in  Figure  B4.  Go  to 
Appendix  C,  Table  Cl  and  select  the  coefficients  a<?,  aj,  and  G.  Enter  these  into 
the  models  in  Figures  B3  and  B4.  Never  use  the  maximum  gain  given  in  Table 
Cl,  but  operate  between  0.01  and  0.1  below  G  or  the  optical  circuit  will  tend  to 
overload. 

c.  The  SIMULINK  subsystem  of  the  photodetector  is  shown  in  Figure  B5.  The 
trigger  frequency  block  is  the  same  as  the  PRF  set  in  the  “Mode-Locked  Laser” 
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block.  The  “Windowed  Integrator”  should  have  an  “Input  vector  size”  of  “1”  and 
an  “Integration  window  length” 


WL  = - - - . 

2  (PRF) 


(Al) 


The  “Sample  time”  is  the  computer  time  steps  and  the  “Method”  of  integration  is 
“First  Order”.  The  “Envelope  Detector”  block  is  set  to  detect  the  “upper” 
boundary  side. 

7.  The  “Output  Comparator”  (Fig.  Bl)  is  a  SIMULINK  relay  function  block.  The 
threshold  should  be  set  at  0.5.  It  should  be  set  to  output  “1”  when  the  block  is  “on” 
and  a  “-1”  when  the  block  is  “off’. 

8.  The  “Computer  Downsample”  SIMULINK  block  (Fig.  B 1)  is  used  to  downsample 
the  computer  time  steps  to  match  the  PRF  prior  to  decimation.  In  each  block,  enter 
the  input  sample  time  and  the  downsampled  integer  factor.  Downsampling  by  factors 
greater  than  8  with  one  function  block  should  be  avoided  due  to  inaccuracies  in  the 
program.  Therefore,  multiple  blocks  are  used. 

9.  The  “Decimation”  SIMULINK  function  block  (Fig.  Bl)  lowpass  filters  and  resamples 
the  signal  to  remove  the  modulation  and  out-of-band  noise.  The  input  sample  time  is 
equal  to  the  PRF.  The  decimation  factor  is  the  required  factor  to  downsample  the 
oversampled  signal  at  the  Nyquist  rate.  As  with  the  downsample  block, 
downsampling  factors  greater  than  8  should  be  avoided.  Again,  use  multiple  blocks 
in  series.  The  lowpass  filter  type  is  “firl”.  Filter  order  of  15-25  is  sufficient  for 
sinusoidal  signals. 

10.  The  “Gain”  function  block  (Fig.  Al)  is  to  increase  the  output  signal  after  filtering  and 
the  scope  is  to  view  the  output  digitized  signal  in  the  time-domain. 
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APPENDIX  C:  PROGRAMS  FOR  COMPARISON  OF  THE  ALL- 
ELECTRONIC  AND  OPTICAL  ACCUMULATORS 

This  appendix  shows  the  programs  used  to  convert  between  the  optical  and  all-electronic 
coefficients  for  comparison  and  contrast  of  the  optical  and  electronic  accumulators. 


%  File :  coef_calc .  m 

%While  loop  to  do  repeated  calculations 

pick=l; 

while  pick=l 

fprintf ( 1  Convert  from  optical  to  electronic  or  vice  versa?\n')  ; 
fprintf  ( '  1-Convert  from  electronic  to  optical  coef fiencients  .  \n T )  ; 
fprint'f  ( r  2-Convert  from  optical  to  electronic  coef  fiencients .  \n' )  ; 
choice=input ( 1  Pick  one:  T); 

%If  know  values  for  electronic  and  want  to  see  how  compares  to  optical 

if  choice==l 

C=input ( 1  Enter  numerator  coefficient  for  H21  or  H12:  1 ); 

D=input ( r Enter  denominator  coefficient  for  H21  or  H12:  ' ); 

CD=C/D; 

fprintf ( 1  Choose  which  value  to  fix\nf); 
fprintf ('1  -  fix  aO\n'); 
fprintf (f2  -  fix  al\nf); 
choose=input ( 1  Pick  one:  '); 

%Must  choose  one  of  the  directional  coupler  coefficients  to  solve 

equation 

if  choose==l 

aO=input ( 1  Choose  a  coefficient  for  aO  (0  —  >1) :  f); 

al=l/ ( (CD*aO/ (1-aO) ) +1) 

G=D/ (a0*al) 

elseif  choose==2 

al=input  (' Choose  a  coefficient  for  al  (0 — >1)  :  ' ); 
a0=l/( (CD*al/(l-al) )+l) 

G=D/(aO*al) 

end 

%If  you  know  optical  coefficients  a.nd  want  to  compare  to  electronic 

coefficients 

elseif  choice==2 

a0=input (' Enter  coefficient  for  a0(0 — >1) :  *); 
al=input ( 'Enter  coefficient  for  al(0 — >1) :  '); 

G=input { 'Enter  Value  for  Gain  (G)  :  '); 

C= ( 1-aO) * (1-al) *G 

D=aO*al*G 

end 

pick=input  ( 1  Type  1  to  repeat  or  2  to  quit:  '); 
end 

Program  Cl :  MATLAB  program  used  to  convert  the  electronic  accumulator  model 
coefficients  to  their  associated  optical  coefficients,,  and  vice  versa. 
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Coupling  Ratio,  aj 


Transfer  Function 


Program  C2:  SIMULINK  program  to  compare  and  contrast  the  optical  fiber  lattice  and 
electronic  accumulators  along  with  their  transfer  function.  Use  program 
Cl  to  calculate  coefficients.  Use  any  input  signal.  (File:  integrator s.mdl) 


Coupling  ratio,  ao 


0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

0.1 

100 

50 

33.33 

25 

20 

16.66 

14.28 

12.5 

11.1 

0.2 

50 

25 

16.65 

12.5 

10 

8.33 

7.15 

6.25 

5.55 

0.3 

33.33 

16.68 

11.11 

8.33 

6.67  '  5.55 

4.76 

4.166 

3.703 

0.4 

25 

12.5 

8.34 

6.25 

5 

4.166 

3.57 

3.125 

2.78 

0.5 

20 

10 

6.68 

5 

4 

3.333 

2.85 

2.5 

2.22 

0.6 

16.66 

8.33 

5.55 

4.17 

3.33 

2.78 

2.38 

2.08 

1.85 

0.7 

14.28 

7.14 

4.76 

3.57 

2.86 

2.38 

2.04 

1.785 

1.587 

0.8 

12.5 

6.25 

4.16 

3.125 

2.5 

2.08 

1.785 

1.563 

1.388 

0.9 

11.1 

5.56 

3.7 

2.78 

2.22 

1.852 

1.588 

1.389 

1.234 

Table  Cl :  Matched  optical  coefficients  for  monotonically  increasing  response  in  the 
fiber  lattice  structure.  By  selecting  a  coupling  ratio  set,  the  correct  optical 
gain  (G)  can  be  found  for  fiber  lattice  monotonic  operation. 
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APPENDIX  D:  PHASE  COHERENT  SECOND-ORDER  OPTICAL 
IA  DIGITAL  ANTENNA 


This  appendix  shows  the  proposed  second  order  EA  digital  antenna 


m 
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APPENDIX  E:  MACH-ZEHNDER  INTERFEROMETER 

MAGNITUDE  AND  DIRECTION  PROGRAMS 

This  appendix  shows  two  separate  programs,  one  MATLAB  (Prog.  El)  and  one 
SIMULINK  (Prog.  E2)  that  model  the  MZI  direction  and  magnitude  intensity  transfer 
functions  (Eqn.  9). 


%  File:  inf plot. m 

G  =  3.0e-6; 
lambda  =  0.9e-6; 
n  =  2.205; 
r  =  30.8e-12; 
gamma  =  1 ; 

L  =  0.0020; 


%inter-electrode  gap  (m) 

%free  space  optical  wavelength  (m) 
%index  of  refraction 
%electro-optic  coefficient  (v/m) 
%overlap  integral 
%length  (m) 


phase  =  pi; 

dphi  =  (2*pi*n^3*r*gamma*L) / (G*lambda) ; 


v  —  -1.5: .04:1.5; 


Imag  =  ( 0 . 5+0 . 5*cos (dphi*v+phase) ) ; 

Idir  =  (0.5+0. 5*cos (dphi*v-phase/2) ) ; 

plot (v,  Imag, 1 o ? ) 
hold 

plot (v, Idir,  ’  +  ’  )  •  • 

axis ( [-1 .5  1.5  0  1] ) 

xlabel ( ’Normalized  Input  Signal  (Volts)1) 
ylabel ( ’Normalized  input  Intensity’) 
%gtext (’ Direction ’ ) 

%gtext ( ’Magnitude ’ ) 
hold 


Program  El :  MATLAB  program  used  to  illustrate  the  MZI  direction  and  magnitude 
transfer  functions. 


Magnitude 

interferometer 


%  File:  inferplot . mdl 

function  I  out=interfer ( v) 


G=3 . Oe-6; 
wl=l . 5e-6; 
n=2 .205; 
r=30 . 8e-12; 
GAMMA=1 . 7 ; 
L=0 . 0020 ; 
phase=pi; 


%  Inter-elect  gap  (m) 
%  Wavelength  (m) 

%  Index  of  refraction 
%  E-Opt  coef  (v/m) 

%  Overlap  integral 
%  length  (m) 

%  Mag=pi,  Dir=~pi/2 


K= (2*pi*nA3*r*GAMMA*L) / (G*wl) ; 
I_out=( .5+.5*cos (K*v  +  phase)); 


Program  E2:  SIMULINK  program  and  associated  MATLAB  function  for  calculating 
the  magnitude  and  direction  MZI  transfer  functions.  The  MATLAB 
function  name  is  inserted  into  the  appropriate  SIMULINK  “MATLAB 
Function”  block. 
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APPENDIX  F:  PHASE  COHERENT  SIMULINK  PROGRAMS  FOR 
MACH-ZEHNDER  INTERFEROMETERS  (MZIs) 

This  appendix  shows  the  SIMULINK  programs  to  compare  and  contrast  the  MZI  phase 
coherent  model  (Eqn.  11)  with  the  MZI  intensity  transfer  function  (Eqn.  4).  This 
program  recreates  Figure  13.  (File:  MZI_phc.mdl) 


Figure  F 1 :  SIMULINK  model  demonstrating  the  phase  coherent  MZI 
versus  the  MZI  intensity  transfer  function. 


Signal  (V)  ri~~1 


Figure  F2:  General  SIMULINK  model  of  a  phase  coherent  MZI. 

Used  for  both  directional  and  magnitude  results.  This 
subsystem  is  part  of  the  MZI  subsystems  in  Figure  FI. 


Figure  F3:  SIMULINK  subsystem  of  the  "phase  modulator"  in 
Figure  F2.  Changing  the  DC  bias  in  this  example  to 
7i  creates  the  magnitude  MZI. 
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Detector 


Figure  F4:  SIMULINK  subsystem  of  phase  coherent  MZI  photodetector  and 
MZI  intensity  function  in  the  MATLAB  function  block.  Either 
output  can  be  chosen  to  go  to  the  scope  (Fig.  FI)  for  viewing. 
This  model  is  used  for  both  the  magnitude  and  direction  MZIs. 
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